NEWTONIAN NORMAL SHIFT 
IN MULTIDIMENSIONAL RIEMANNIAN GEOMETRY. 



RUSLAN A. SHARIPOV 

Abstract. Explicit description for arbitrary Newtonian dynamical system admit- 
ting the normal shift in Riemannian manifold of the dimension n ^ 3 is found. On the 
base of this result the kinematics of normal shift of hypersurfaces along trajectories 
of such system is studied. 



1. Introduction. 

In series of papers [1-16] written in 1993-1996 a theory was constructed that 
determines and describes special class of Newtonian dynamical systems ad- 
mitting the normal shift of hypersurfaces in Riemannian and Finslerian mani- 
folds. On the base of these papers two theses were prepared: thesis for the degree 
of Doctor of Sciences in Russia [17] and thesis for the degree of Candidate 
of Sciences in Russia [18]. However, some results included in thesis [17] are still 
not published in journals (see [10] and [11]). 

Moreover, when preparing thesis [17], in paper [16] an error was found. Eliminat- 
ing this error led to new result that consists in complete and exhausting description 
of all Newtonian dynamical systems admitting the normal shift in Riemannian 1 
manifolds of the dimension n ^ 3. The goal of this paper is to explain this new 
result, and to give a description for kinematics of normal shift of hypersurfaces, 
more detailed than it was possible before now. On the base of the same result one 
can get new (more simple) proof for the main theorem from unpublished paper 
[11] (see also § 7 in Chapter VI of thesis [17]), and one can answer the question by 
A. V. Bolsinov and A. T. Fomenko, which they asked when author reported thesis 
[17] in the seminar of the Chair of Differential Geometry and its Applications at 
Moscow State University. 

Classical construction of normal shift of hypersurfaces in Riemannian manifold 
M is well-known. In its original form it arises in the case when M is R 3 with 
standard flat Euclidean metric. Let S be two-dimensional surface in K 3 . From 
each point p on S we draw a segment of straight line in the direction of normal 
vector n = n(p). Denote by p t the second end of this segment. When p runs over 
5 1 , point pt sweeps some other surface St as shown on Fig. 1.1 below. So we have 
the map ft'. S — > St known as classical normal shift or as Bonnet transformation. 

When transferring from M = R 3 to the case of arbitrary Riemannian manifold 
M, we replace surfaces by hypersurfaces, and straight line segments, connecting 

x This result has no direct generalization for the case of dynamical systems in Finslerian ma- 
nifolds (this case was considered in [10], see also Chapter VIII in thesis [17]). But, nevertheless, 
it possibly has some analog in Finslerian geometry. This problem is not yet studied. 
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p and pt, by the segments of geodesic lines. In this form classical construction of 
normal shift of hypersurfaces is known as geodesic normal shift. Construction of 

geodesic normal shift contains a numeric pa- 
rameter t, i. e. it determines the whole family 
of hypersurfaces St- When parameter t is var- 
ied, point p t moves along geodesic lines (here 
they are called trajectories of the shift). In 
local coordinates x 1 , . . . , x n on M they are 
described by ordinary differential equations 

n n 

(i.i) r • XX't 

i=l j=l 

where k = 1, . . . , n. By r*. = -.(x 1 , . . . ,x n ) 
in (1.1) we denote components of standard metric connection V for Riemannian 
metric g on M. The property of normality of geodesic shift is expressed by the 
following well-known fact. 

Theorem 1.1. All hypersurfaces St in the construction of geodesic normal shift 
are perpendicular to trajectories of shift. 

In other words, trajectories of shift described by the equations (1.1) cross each 
hypersurface S t transversally; at the points of intersection they pass along normal 
vectors to St- 

The idea of generalizing the construction of normal shift, which was realized 
in papers [1-16], is very simple. It consists in replacing (1.1) by slightly more 
complicated ordinary differential equations in t: 

n n 

(i.2) ■ Y^l v -j y}J F V *")■ 

»=i j=i 

When M = R 3 and = 0, the equations (1.2) express Newton's second law: they 
describe the motion of a mass point with unit mass in the force field determined by 
right hand sides of these equations. In the case of arbitrary Riemannian manifold M 
these equations, as appears, also have physical interpretation. They describe the 
dynamics of complex mechanical systems with holonomic constraints. Manifold 
M arises as configuration space of such systems, its dimension is determined by 
actual number of degrees of freedom (upon resolving all constraints). Thereby M 
is canonically equipped with the structure of Riemannian manifold, its metric is 
given by quadratic form of kinetic energy: 

^ n n 
i=l j=l 

(see details in Chapter II of thesis [17]). Due to the analogy with Newton's second 
law the equations (1.2) arc called the equations of Newtonian dynamical system 
on Riemannian manifold. Vector F, whose components are given by right hand 
sides of the equations (1.2), is called a force vector. It determines force field of 
Newtonian dynamical system (1.2). 
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Note that the choice of local coordinates in defining Newtonian dynamical system 
(1.2) is of no importance. By the change of local coordinates the shape of the equa- 
tions remains unchanged, thought the components of connections and components 
of force vector are transformed according to standard formulas, which are well- 
known from course of differential geometry (see [19-22]). This property expresses 
coordinate covariance of differential equations (1.2). 

2. Newtonian normal shift of hypersurfaces. 

Having formulated the idea of generalizing the construction of geodesic normal 
shift, we shall describe how it was realized in papers [1-16]. Let S be some hyper- 
surface in M and let p be some point on S. In local coordinates x 1 , ... , x n on M 
such point p is characterized by its coordinates a; 1 ^), ... , x n (p) and by normal 
vector n(p) at this point. Let's use Newtonian dynamical system (1.2) in order to 
define a shift of hypersurface S. With this aim let's associate each point p G S with 
the following initial data for the system of differential equations (1.2): 

(2.1) x k \ = x k {p), x k \ =v{p)-n k {p). 

\t=o lt=o 

Here n k {p) is fc-th component of normal vector n(p), and v{p) is some scalar quan- 
tity depending on the point p G S. Solving 
Cauchy problem with initial data (2.1) for the 
equations (1.2), we obtain a set of functions 



(2.2) 



' x 1 = x 1 ^^), 



k x n = x n (t,p). 



These functions define in parametric form the 
trajectory r = r(t,p) of Newtonian dynamical 
system with force field F. This trajectory at 
initial instant of time t = crosses hypersur- 
face S at the point p, passing in the direction 
of unitary normal vector n(p). Parameter v(p) 
in initial data (2.1) determines the modulus of initial velocity for this trajectory: 

(2.3) v| =i/(p).n(p). 

\t=o 

Choice of local coordinates in defining trajectory r = r(t,p) is of no importance. 
Change of local coordinates changes functions (2.2), but it doesn't change the curve 
r = r(t,p). This is due to coordinate covariance of differential equations (1.2) and 
coordinate covariance of initial data (2.1). 

Drawing trajectories r = r(t,p) outgoing from all points peS and taking points 
p t = r(t,p) that corresponds to some fixed value of parameter t, we obtain the 
hypersurface St and displacement map f t : S — > St- However, we should remember 
two nuances. Parameter t for trajectory r — r(t,p) of dynamical system (1.2) do 
not coincide with its length. The range of this parameter t always includes initial 
point t = 0, but it can be a restricted interval 



hip) < t < t 2 [p). 
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Upper and lower bounds of this interval in general case depend on the point p G S. 
Hence for a fixed value of t the displacement map f t : S — > St can be defined not 
for all points p G S. 

Second nuance is due to singular points (caustics) that may appear on the hyper- 
surface S t for large enough values of parameter t. This imposes one more restriction 
onto the range of parameter t. Note that this restriction is present in classical con- 
struction of geodesic normal shift as well. 

The above two nuances restrict possible range of parameter t. However, if we 
arc interested in small values of t only (as below) , we can use the following lemma. 

Lemma 2.1. If parameter v(p) in (2.3) is a smooth nonzero function on the hy- 
persurface S, then for each p £ S there exists some neighborhood S' — Os{p) on S 
and there exists a number e such that displacement maps ft : S' — > S^ are defined 
for all t G (— e, +s). They form smooth one-parametric family of diffeomorphisms. 

This lemma is an immediate consequence of theorem on existence, uniqueness, 
and smooth dependence of initial data for the solution of Cauchy problem for sys- 
tems of ODE's (see [23] and [24]). Taking into account lemma 2.1, we can consider 
displacement maps ft : S — > St, which are possibly defined only locally on S, as 
a construction of shift of hypersurface S along trajectories of Newtonian dy- 
namical system with force field F. Function v(p) on S is a parameter in such 
construction of shift. 

Shift of hypersurface S along trajectories of dynamical system (1.2), as it was 
constructed above, possess the property of normality at the initial instant of time 
t = 0. This means that trajectories of shift are passing through initial hypersurface 
S along normal vectors on it. Does this property persist for t ^ 0, i. e. can we 
prove theorem similar to theorem 1.1 ? The answer to this question in general case 
is negative (see examples in [18]). But there are special cases, when the property 
of normality persists for all instants of time. We describe them by formulating the 
following definition. 

Definition 2.1. Shift ft - S — > St of hypersurface S along trajectories of Newtonian 
dynamical system with force field F is called a normal shift if all hypersurfaces 
St (for all permissible values of parameter t) are orthogonal to trajectories of shift. 

3. Dynamical systems admitting the normal shift. 

What does the property of normality for shift f t : S — ► St depend on ? On the 
choice of hypersurface S ? On the force field F of Newtonian dynamical system ? 
We also have the opportunity to choose the function v(p) on S. In the case of 
identically zero force field F = (which corresponds to geodesic flow on M) the 
choice v(p) = 1 provides normality condition from definition 2.1 for arbitrary initial 
hypersurface S C M. Are there some other force fields with similar property? The 
aim to know this was the motivation for writing preprint [1]. In this preprint 
we introduced the concept of Newtonian dynamical system admitting the 
normal shift. This concept has become a central point for later investigations. 

Definition 3.1. Newtonian dynamical system on Ricmannian manifold M is called 
a system admitting the normal shift if for any hypersurface S in M, and for 

any point p G S, there is a neighborhood S' — Os(po) of the point p n on S, and 
there is a smooth function v{p) in S', such that the shift f t : S' — > S' t defined by 
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the function v(p) is a normal shift along trajectories of considered dynamical in the 
sense of definition 2.1. 

The condition stated in definition 3.1 was called the normality condition for 

Newtonian dynamical system with force field F. First we considered the case M — 
M 2 (see [1]). In [1] we derived partial differential equations for the components of 
force field F which, when being fulfilled, are sufficient to provide normality condition 
from definition 3.1. These equations were called the normality equations or, more 
exactly, weak normality equations. In preprint [1] we also constructed first non- 
trivial examples of dynamical systems that admit normal shift. When generalizing 
these results from M = W 2 to multidimensional case M = R n in preprint [1], we 
have found that weak normality equations should be supplemented by so called 
additional normality equations. All above results in brief form were announced in 
[4] . Their full version were published in [2] and [3] . Later in papers [6] and [7] they 
were generalized for the case of arbitrary Riemannian manifold M . Main result of 
paper [6] is the derivation of weak normality equations for this more complicated 
geometric situation. We write these equations without comments so far: 



(3.1) 



3 = 1 



E E ( v ^- + v ^ - 2 v ^ F * F o) Ni n + 

»=i j=i 



i=l j=l \ r =l / 



Additional normality equations for the force field F of Newtonian dynamical sys- 
tems on Riemannian manifolds were derived in [7] . They look like 



(3.2) 



run 



i=i j=i 



k m— 1 



EE^^fE^^r^-v^, 



i=l j=l 



s m— 1 



n n n 



EE'^,'"'' EEE i 

i—l j—1 i= 1 j—1 m=l 



n-l 



Results of papers [6] and [7] were announced in [12]. In deriving the equations (3.2) 
we have found that it is convenient to make slight modification of definition 3.2. 
We added normalizing condition for the function v(jp): 



(3.3) 



v(p ) = vq. 



As a result we obtained strong normality condition. It is formulated as follows. 
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Definition 3.2. Newtonian dynamical system on Ricmannian manifold M is called 
a system admitting the normal shift in strong sense if for any hypersurface S in 
M, for any point po G S, and for any real number v ^0 there is a neighborhood 
S' = Os(po) of the point po on S, and there is a smooth function v{p) in S' 
normalized by the condition (3.3), such that the shift f t : S' — ► S' t defined by the 
function v{p) is a normal shift along trajectories of considered dynamical in the 
sense of the definition 2.1. 

Strong normality condition from definition 3.2 implies normality condition for- 
mulated in definition 3.1. The relation of strong normality condition with normality 
equations (3.1) and (3.2) is described by the following theorem. 

Theorem 3.1. Newtonian dynamical system on Riemannian manifold M admits 
the normal shift in strong sense if and only if its force field F satisfies normality 
equations (3.1) and (3.2) simultaneously. 

Theorem 3.1 was proved in [13]. Detailed version of this proof can be found in 
thesis [17]. We shall not give this proof here, since this would require to reproduce 
many details of derivation of the normality equations (3.1) and (3.2), and hence 
this would be doubling for the papers [6] and [7]. Instead, we shall give comments 
to normality equations (3.1) and (3.2). 

4. Extended algebra of tensor fields. 

Normality equations (3.1) and (3.2), as well as the equations (1.2), possess the 
property of coordinate covariance. Writing these equations, we assume that 
some local coordinates x 1 , . . . , x n in M are chosen. Under the change of local 
coordinates all quantities, which are contained in the equations (3.1) and (3.2), 
do change according some definite rules. However, this do not change the shape 
of these equations. Components of force vector F in (3.1) and (3.2) depend on 
x 1 , ... , x n , and on components of velocity vector v; the latter is a tangent vector 
at the point p with coordinates x 1 , ... , x n . This means that vector F depend on 
the point q — (p, v) of the tangent bundle TM. 

Definition 4.1. Vector- function F that for each point q = (p, v) of tangent bundle 
TM puts into the correspondence some vector from tangent space T P (M) at the 
point p = n(q) is called an extended vector field on the manifold M. 

Here tt : TM — ► M is a map of canonical projection from TM to the base 
manifold M. Let's consider the following tensor product: 

r times 

TJCp, M ) = T P (M)®~.®T P (M)® T p *(M) ® . . . ® T p *(M) . 

V * ' 

s times 

Linear space XJ(p, M) is called a space of tensors of the type (r, s) at the point 
p of the manifold M . Elements of this space are called r-times contravariant and 
s-times covariant tensors, or simply tensors of type (r, s) at the point p € M. 

Definition 4.2. Tensor-valued function X that for each point q of tangent bundle 
TM puts into the correspondence some tensor from the space TJ(p, M) at the point 
p = w(q) is called an extended tensor field of type (r, s) on the manifold M. 



NEWTONIAN NORMAL SHIFT 



7 



In local coordinates x 1 , ... , x n on the manifold M extended tensor fields are 
expressed by the functions of double-set of arguments: 

(4.1) Xl-;X(x\... ,x n ,v\... ,v n ). 

In the normality equations (3.1) and (3.2) we can see components of several ex- 
tended tensor fields. Velocity vector v by itself can be considered as extended 
vector field on M. Its modulus v = |v| is an extended scalar. Extended vector field 
N with components TV 1 , . . . , N n , which are contained in the equations (3.1) and 

(3.2) , is defined as the following quotient: 

v v 

(4.2) N=- = — . 

v |v| 

This is the field of unitary vectors directed along the vector of velocity. And finally, 
in equations (3.1) and (3.2) we have components of operator field P. This is the 
field of orthogonal projectors onto the hyperplane perpendicular to the velocity 
vector. Components Pj of this field are given by the formula 

(4.3) P i j =5)-N i N j . 

Let TJ(M) be the set of smooth extended tensor fields of type (r, s). This set 
has the structure of module over the ring of extended scalar fields. The following 
sum is a graded algebra over this ring with respect to tensorial multiplication: 

oc oc 

(4.4) t(m) = 00t;(m). 

Algebra (4.4) is called an extended algebra of tensor fields on the manifold 
M. In extended algebra of tensor fields T(M) one can define two operations of 
covariant differentiation, we denote them by V and V: 

v : t;(m) - r; +1 (M), v : t;(m) - t; +1 (m). 

In local coordinates the result of applying covariant differentiation V to a tensor 
field X with components (4.1) is expressed by the following formula: 

nyh...i r n n ayii-.-ir 

v Y H...i r _ uyv ji-j, \^\^„a r b h-3s , 



(4.5) 



dx m ma dv b 

a=l 6=1 

r n s n 

i \ ^ pifc yii--- a k ... i r yb k yii 

~ Z—i Z^i ma k ji js L^i mjk j 



mj k j 1 ...b k ...j s 

k=l a k = l fc=l b fc = 1 



The result of applying V to X is expressed by less complicated formula: 
(4-6) V m A-- = —i-^. 

Formula (4.6) for the components of the field VX contains only derivatives with 
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respect to components of velocity vector. Therefore V is called velocity gradient. 
Covariant differentiation V defined by formula (4.5) is called spatial gradient. 

Defining operators V and V by means of formulas (4.5) and (4.6), we assume 
that some local coordinates are chosen. This way of defining V and V is quite 
sufficient for our purposes in the theory of newtonian dynamical systems admitting 
the normal shift. But there is another (invariant) way of defining these operators, 
ft is based on the analysis of differentiations in the extended algebra of tensor fields. 

Definition 4.3. The map D : T(M) -> T(M) is called a differentiation in 

extended algebra of tensor fields if the following conditions are fulfilled: 

(f) compatibility with grading: £>(TJ(M)) C 7J(M); 

(2) R-lincarity: D(X + Y) = D(X) + £>(Y) and D(AX) = AL>(X) for A e R; 

(3) pcrmutability with contractions: D(C(X)) = C(D(X.)); 

(4) Leibniz rule: £>(X <g> Y) = D(X) <8> Y + X <8> £>(Y). 

Among the results of thesis [17] the following structural theorem is worth to 
mention here. It describes the structure of all differentiations in extended algebra 
of tensor fields T(M). 

Theorem 4.1. Let M be smooth real manifold equipped with some extended affine 
connection T . Then each differentiation D in extended algebra of tensor fields T(M) 
on this manifold breaks into the sum 



where X and Y are some extended vector fields, and S is a degenerate differentiation 
defined by some extended tensor field S of type (1, 1) in M. 

Theorem 4.1 is an analog of structural theorem for differentiations in the algebra 
of ordinary (not extended) tensor fields (see [19]). 

5. Reduction of normality equations in the dimension n ^ 3. 

If we take into account formulas (4.2), (4.3), (4.5), and (4.6), we see that nor- 
mality equations (3.1) and (3.2) form strongly overdetermined system of partial 
differential equations with respect to components of force field F of Newtonian 
dynamical system. Analysis of this system (see [16]), is based on scalar ansatz 



Formula (5.1) expresses components of force vector through one scalar field A, 
which is interpreted as the projection of F onto the direction of velocity vector. 
This formula follows from first part of equations in the system (3.1). Therefore, 
when substituting (5.1) into weak normality equations (3.1), first part of these 
equations appears to be identically fulfilled. While second part is brought to 



£ , = Vx + V Y + S 



n 



(5.1) 




(5.2) 




= 0. 
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For n = 2 the equations (5.2) exhaust whole list of reduced normality equations. 
The matter is that in two-dimensional case, as we mentioned above, additional 
normality equations do not arise at all. While the equations (5.2) are reduced to 
the only one nonlinear partial differential equation for the function A(x 1 ,x 2 ,v 1 , v 2 ). 
Detailed study of this equation is given in the thesis by A. Yu. Boldin [18]. 

In multidimensional case n > 3 the process of reducing normality equations can 
be moved much further. Substituting (5.1) into the first part of additional normality 
equations brings them to the following form: 

n n / n \ 

Yl P ° P e V ^ S A + Y V r A W V q V s A = 

(5.3) 8=1 r=1 V 9=1 / 

= E E P ° P e V °^ r A + Y V*A N q V q V r A . 

s=l r=l \ q=l / 

Similarly, substituting (5.1) into second part of the equations (3.2) gives 

n n 

(5.4) ^^F;v r V s AF £ = A^. 

r— 1 s— 1 

Here A is a scalar quantity, the value of which is uniquely determined by the equa- 
tions (5.4) even if we do not know it a priori: 

. ff P rS V r V s A 

x -1^2^ n _ 1 ■ 

r—l s—1 

The equations (5.4) are most remarkable. According to the formula (4.6) they con- 
tain only the derivatives with respect to the variables v 1 , ... , v n . This corresponds 
to varying the function A(x , ... , x n , v , . . . , v n ) within the fiber of TM over the 
fixed point with coordinates x , . . . , x n in the base manifold M. 

Definition 5.1. Extended tensor field X on Riemannian manifold M is called 
fiberwise spherically symmetric if |vi| = |v 2 | implies X(p, vi) = X(p, v 2 ). 

In other words, fiberwise spherically symmetric extended tensor fields depend 
only on modulus of velocity vector within fibers of tangent bundle TM. Such fields 
naturally arises in the analysis of the equations (5.4). Here we have the following 
theorem proved in paper [16]. 

Theorem 5.1. Extended Scalar field A on Riemannian manifold M satisfies equa- 
tions (5.4) if and only if it is given by formula 

n 

(5.5) A = a + y^bjV*, 

where a is some fiberwise spherically symmetric scalar field, and bi are components 
of some fiberwise spherically symmetric covectorial field b. 

Further substitution of (5.5) into the equations (5.2) and (5.3) yields the following 
equations with respect to fields a and b: 
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Here a = a(x 1 , ... , x n ,v) and bi = b^x 1 , . . . ,x n ,v). Variable v denotes the modu- 
lus of velocity vector: v = |v|. 

Equations (5.6) and (5.7) is bound with most dramatic instant in constructing 
theory of dynamical systems admitting the normal shift. By deriving these equa- 
tions in paper [16] the mistake was made, nonlinear terms in (5.7) were omitted. 
As a result (5.7) looked like db r /dx s = db s /dx r . Further analysis of erroneous 
equations has led to the ordinary differential equation 

(5.8) y" = H y (y' + 1) + H x 

where H = H(x,y), H y = dH/dy, H x = dH/dx. With the aim to find as more 
functions H(x,y), for which the equation (5.8) is explicitly solvable, as possible we 
considered the following change of variables: 

(5.9) = 

I y = y{x,y). 

The equations that could be brought to the form (5.8) by means of change of 
variables (5.9) belong to the following class of equations: 

(5.10) y" = P(x, y) + 3 Q(x, y)y' + 3 R(x, y) y' 2 + S(x, y) y' 3 . 

Study of point transformations in the class of equations (5.10) has the long his- 
tory (see [25-46]). However, we couldn't find an answer to the question: how to 
extract the equations (5.10) that could be brought to the form (5.8) by means of 
point transformation (5.9). This stimulated our own investigations (see [47-50]). 
We managed to get some results in describing classes of point equivalence for the 
equations (5.10). But now, since the error in [16] is found, these results have sepa- 
rate value, which is not related to the theory of dynamical systems admitting the 
normal shift. And we are to return to the equations (5.6) and (5.7). 

6. Derivation of reduced normality equations. 

For the beginning let's derive the (5.6) and (5.7) by substituting (5.5) into the 
equations (5.2) and (5.3). Denote by a' and 6- the derivatives 

a' - — b '■ - — 

dv 1 1 dv 

Let's do the calculations necessary for substituting (5.5) into (5.2) and (5.3): 

n 

(6.1) V S A = V.al^V.M', 

i=l 

(6.2) V r V s ^ = (v r a' + J2 ^rK v^j N s + V r b s . 



NEWTONIAN NORMAL SHIFT ... 11 



(6.3) V a A=^a' + J2V i v i jN a + b a , 

V r V s A = ( a" + b" v l J N r N s + 

(6.4) V i=i / 

+ b' s N r + b' r N s + ( - + K Ni J P rs- 

From formulas (6.3) and (6.4) for derivatives we obtain the following relations: 

n n 

(6.5) £p;v r A = £p;6 r , 

r—1 r—1 

n n n 

(6.6) E P e m = E P e K- 

s—1 q—1 s—1 

Let's combine (6.5) and (6.6). As a result we get the relationship 

n n n n n 

EEE P ^ P ° N q V q V s A = E K P e b r K- 

s—1 r—1 q—1 s—1 r—1 

Then let's multiply both sides of the relationship (6.2) by P e s and contract with 
respect to pair of indices r and s. This yields one more relationship: 

n n n n 

]T ]T PI PI V r \7 s A = E P « P e Vrh. 

s — 1 r—1 s — 1 r—1 

Now, if we add above two relationships, we get the result of substituting (5.5) into 
the left hand side of the equation (5.3): 

n n / n \ 

E p * p e v r v s A + v r A N q v q v s A = 

fn rj\ 8=1 r=l \ q=l / 

s—1 r—1 s—1 r—1 

Similarly we calculate right hand side of the equation (5.3): 

n n / n \ 

E P ° P e ( ^sV r A + V S A N q \7 q V r A ) = 

(6.8) 



s — 1 r—1 \ 9—1 

n n n n 



s—1 r—1 s—1 r—1 

On the base of (6.7) and (6.8) we conclude that the equation (5.3) is reduced to 

n n 

( 6 - 9 ) E E P ° P e i^rbs + b r b' s - V s b r - b s b' r ) = 0. 

3=1 r—1 
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For the further analysis of the obtained equations (6.9) one should use the peculiar- 
ity of covectorial field b from extended algebra of tensor fields on M. Components 
of this field b\, . . . , b n depend only on modulus of velocity vector, but they do not 
depend on its direction. By calculating derivatives V r & s and V s b r in (6.9) we apply 
the following theorem. 

Theorem 6.1. Let X] 1 •"*. r = X] 1 '" l i r (x 1 , . . . ,x n ,v) be components of fiberwise 

J 1 • • • j S J 1 ' ' ' J s 

spherically symmetric tensor field X from extended algebra T(M). Then compo- 
nents of spatial gradient VX for this field are given by formula 

Yj yii — ir _ h---3s i \ 1 \ 1 -pit V* 1 " - *r _ 
Vm ^jl-js - Q x m ' ..' 3s 

(6.10) ' ' 



s n 

E\ ^ I vh 
Z^ mj k ji...b k . 
k=l b k = l 



Proof. Formula (6.10) is obtained as a result of reduction from formula (4.5). For 
the components of fiberwise spherically symmetric tensor field X (as in the state- 
ment of theorem) natural arguments and v, where 



(6.11) u=|v| = 



n n 



1 j=l 



While partial derivatives in formula (4.5) are assumed to be respective to the vari- 
ables x 1 , ... , x n , v 1 , ... , v n . Recalculation of these derivatives to natural vari- 
ables for spherically symmetric field consists in the following substitutions: 

(612) ^S by ^S.^ 

oxi'-i' ox i'-)' oxi'-y dv 
(6-13) by ~^r ± + 



dx m dx m dv dx m ' 

Derivatives dv/dv b and dv/dx m are calculated due to (6.11). Upon finding explicit 
expressions for these derivatives and upon making substitutions (6.12) and (6.13) 
in formula (4.5), we get two extra summands: 



^ ^ dv 2 dx m v ^ ^ ma dv 

a=l b=l a=l b=l 

If we take into account the explicit formula for components of metric connection 
f fi -m\ r b _ 1 n b Q ( I ^Smg dg ma \ 

1 ' ma 2 2_^y \dx m dx a dxl J 

3=1 



(see [19-22]), then we easily see that above two summands cancel each other. As a 
result formula (4.5) transforms into the form (6.10). □ 
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Corollary. Components u rs = V r b s +b r b' s — V s b r — b s b' r of skew- symmetric ex- 
tended tensor field uj in the equations (6.9) depend only on modulus of velocity 
vector |v|, but they do not depend on the direction of the vector v. 

This fact immediately follows from the formula (6.10). It allows us to make 
further simplifications in the equations (6.9). Let c and d be arbitrary two vectors 
from tangent space T p (M). In multidimensional case n ^ 3 we can rotate veloc- 
ity vector v, keeping its modulus unchanged, and can direct it so that it will be 
perpendicular to vectors c and d simultaneously. Then 



a=l e=l 

Therefore the equations (6.9) are transformed as follows: 

n n 

(6.15) E ^ ^ ( Wrbs + br b ' s ~ Vsbr ~ bs ^ = °' 
s— 1 r— 1 

Since c and d are arbitrary two vectors, we can further simplify the obtained 
equations (6.15), bringing them to the form 

(6.16) V r b s + b r b' s = W s b r + b s b' r . 

Next step consists in reducing the equations (5.2). In order to do it we substitute 
(5.5) into (5.2). From (6.1) we derive 

n n n n 

(6.17) Yl V * A P fe = E V * a P k + E E Vsbr v r P s k . 

S — l S—l S—l r — 1 

Then from (6.3) and (6.4) we obtain the following two relationships: 

n n 

(6.18) j2 pqr Vi A = Y. b ° psr > 

r—1 s—l 

n n / n \ 

(6.19) \v\Y l Y l pqr VrV s AP£= {a' + YKv r \P q k . 

r—1 s—l \ r—1 / 

In (6.18) we have free index r, and in (6.19) we have free index q. Let's multiply 
these two equalities (6.18) and (6.19) and do contract with respect to indices r and 
q upon multiplying the resulting equality by g rq . This yields 

n n n / n \ n 

M E E E ^ V^sA PZ=la' + Y K v r W b s P s k . 

q—1 r—1 s—l \ r—1 / s—l 

One more relationship is obtained from (6.4) upon multiplying by N r AP* and 
upon contracting with respect to r and s: 

n n / n \ n 
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Let's multiply (6.2) by N r and P%, then contract it with respect to r and s: 

n n n n 

i v i E E ^ r v ^ = - E E ^ r v ^ p fe • 

r— 1 s — 1 r— 1 s— 1 

Now, in order to write the result of substituting (5.5) into the equations (5.2), we 
have to add (6.17) and three above equalities: 

n / n \ 

V s a + b s a' - a b' s + £ v r (\7 s b r + b s b' r - V r b s - b r b' s ) P fe s = 0. 

s=l \ r=l / 

Let's take into account (6.16), this leads to vanishing the whole expression under 
summation with respect to r. As a result we obtain the following equation: 

n 

(6.20) ^2(V s a + b s a' -ab' s )P£ = 0. 

s=l 

Equations (6.20) are analogous to the equations (6.9), the operation of contraction 
with components of projector P can be omitted: 

(6.21) V s a + b s a' = ab' s . 

Arguments used in deriving the equations (6.21) are similar to those used in deriving 
(6.16) from (6.9). 

Theorem 6.2. Force field F given by scalar ansatz (5.1) corresponds to some New- 
tonian dynamical system admitting the normal shift on the Riemannian manifold 
M if and only if scalar field A in ansatz (5.1) is defined by formula (5.5), while 
extended fields a and b in (5.5) are fiberwise spherically symmetric and satisfying 
the equations (6.16) and (6.21). 

Note that the equations (6.16) and (6.21) coincide with reduced normality equa- 
tions (5.7) and (5.6) we were to derive. 

7. Analysis of reduced equations. 

With the aim of further study of the equations (5.6) and (5.7) let's express 
covariant derivatives in them through partial derivatives. In order to do it we use 
formula (6.10) and take into account symmetry of connection components: 

The equations (7.2) form closed system of equations with respect to components 
b\, ... , b n of covector field b. We can study them separately. Let's consider the 
differential operators in these equations: 

d d 

(7.3) Lj = — t - + b t — , where i=l, ... ,n. 



-(• 


8_\ 

dv J 


)b s , 


/ d 


+ b r 




K dx r 


dv J 
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Now by means of direct calculations we can check that the equations (7.2) are 
exactly the conditions of permutability of operators (7.3): 

(7.4) [L s , L r ] = 0. 

Let R + = (0, +oo) be positive semiaxis on real axis M. Operators (4.3) have 
natural interpretation as vector fields on the direct product of manifolds M x M+. 
Let's complement Li, ... , L„ by one more vector field L„+i, which possibly is 
not commutating with fields Li, ... , L„, but which should complete Li, ... , L„ 
up to a moving frame on the manifold M x K + . Each field Lj has its own local 
one-parametric group of local diffeomorphisms (see [19]) with parameter t^: 

(7.5) tfiiy') : M x R + — >• M x R + . 

Let's fix some point p E M x R + and let's consider composition of such diffeomor- 
phisms applied to the point p : 

(7-6) piy 1 ,... ,y n ,w) = (p 1 (y 1 )° ... °tp n (y n )otp n+1 (w)(p ). 

In left hand side of the equality (7.6) we have the point p parameterized by real 
numbers y 1 , ... , y n , w. This is equivalent to defining local coordinates on M x R + 
in some neighborhood of the point po- Permutability of vector fields (7.4) implies 
permutability of first n maps (7.5) in the composition (7.6). For the vector fields 
Li, . . . , L„ this fact yields the following expressions: 



(7.7) 



d ^ dx % d dv d 
dy 1 4^ dy 1 dx % dy 1 dv 1 



d dx l d dv d 

dy n — ^ dy n dx % dy n dv 



Let's compare formulas (7.7) and (7.3) for vector fields Li, ... , L„. This yields 



(7-8) 753 = 4 = 



dx l _ i _ f 1 for i = k, 
'dy T ~ k ~ \ for k. 



From the same comparison for the functions bu in variables y 1 , ... , f, w we get 



(7-9) b k = 



dv 
dy k 



The relationships (7.8) show that newly constructed local coordinates y 1 , ... , y n , w 
on M x M. + and initial local coordinates x 1 , . . . , x n , v on this manifold are related 
by means of only one function V(y 1 , ... ,y n ,w): 



(7.10) 



x 1 = y 1 , ■ ■ ■ , x n = y n , 
v = Viy 1 , . . . ,y n ,w). 
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Inverse relation is also given by the only one function W{x 1 , ... ,x n ,v): 
(7.11) 



y 1 = x 1 , ... , y n = x r ' 



I w = W(x 1 , . . . , x n ,v). 

Functions V(y 1 , . . . ,y n ,w) and W(x 1 1 . . . ,x n ,v) in (7.10) and in (7.11) are bound 
by the obvious relationships that express the fact that the changes of variables 
(7.10) and (7.11) are inverse to each other: 

Vix 1 ,... ,x n ,W(x 1 ,... ,x n ,v)) =v, 

(7.12) 

W(x\... ,x n ,V(x 1 ,... ,x n ,w))=w. 

Let's use the following natural notations for partial derivatives of the first order: 

! „ dV{x 1 ,... ,x n ,w) 
Vi(x ,...,x ,w) = —. , 

(7.13) 

v w ( x \..., x »M = mx '- dw ' xn ' w) . 

Analogous notations will be used for partial derivatives of the second order: 

x d 2 V(x 1 ,...,x n ,w) 

V^X,...^ ,«,) = , 

(7.14) 

V (x 1 x- w) ^Vjx\...^ lW ) 
v lw [x ,...,x ,w) QxiQw 

Now we can rewrite formula (7.9) in initial local coordinates x 1 , . . . , x n , v. 

(7.15) b k = V k (x 1 ,...,x n ,W(x 1 ,...,x n ,v)). 

Theorem 7.1. Functions b\, . . . , b n satisfy nonlinear system of partial differential 
equations (7.2) if and only if they are determined by some function V(x 1 , ... , x n , w) 
with non-zero derivative dV/dw according to the formula (7.15). 

Proof. In theorem 7.1 we have two propositions. Direct proposition is already 
proved: each solution of the system of equations (7.2) is given by (7.15). Conversely, 
suppose that some function V(x 1 ,... ,x n ,w) with non-zero derivative dV/dw is 
chosen. From dV/dw ^ 0, relying on the theorem on implicit functions (see [51], 
[52]), we derive the existence of the function W(x 1 , . . . , x™, v) such that it is bound 
with V(x 1 7 ... ,x n ,w) by the relationships (7.12). Differentiating these relation- 
ships and taking into account the notations (7.13), we derive 

(7.16) 3W 



(7.17) 



dv V w (x x ,... ^"sW^x 1 ,... ,x n ,v)) 

dW _ V s (x\... ,x n ,W(x\... ,x n ,v)) 
Sx 7 " Vwix 1 , . . . ,x™, iy(x\ . . . ,x n ,v)) 
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Let's substitute V and W into (7.15) and calculate functions bi, ... , b n . Then by 
means of direct calculations we check that the functions obtained satisfy differential 
equations (7.2). Indeed, here we have 



(7.18) 



(7.19) 



^ = V rs (x\... ,x n ,W(x\... ,x n ,v)) + 

+ V rw (x\...,x n ,W(x 1 ,...,x n ,v))—, 

dh 

b s -^ = V s (x\... ,x n ,W(x\... ,x n ,v))x 

dW 

x V rw (x\... ,x n ,W(x\... ,x n ,v))—. 



Let's add the equalities (7.18) and (7.19) and let's take into account formulas (7.16) 
and (7.17) for derivatives. This yields the equality 

(7.20) (J^+b s -^b r = V rs (x 1 ,... ,x n ,W(x\... ,x n ,v)). 

Due to first formula (7.14) we can transpose indices r and s in right hand side of 
(7.20), i. e. V rs = V sr . This provides the equations (7.2) for the functions (7.15) we 
have constructed above. □ 

Having constructed general solution for the equations (7.2), now let's study the 
equations (7.1). Let's implement the change of variables (7.10) and transfer to 
variables y 1 , . . . , y n , w. In left hand side of (7.1) we have the same differential 
operator L s as in (7.2). In variables y 1 , ... , y n , w this operator is written as 
L s = d/dy s (see relationships (7.7)). Let's transform the operator in right hand 
side of (7.1) to the variables y 1 , . . . , y n , w: 

d ^ dy l d dw d dW d 

dv ^— ' dv dv 1 dv dw dv dw 

»=i y 

For further transformation of the above expression for operator d/dvwe use formula 
(7.16). As a result we get the following relationship: 

d 1 d 



dv Vwiy 1 , ... ,y n ,w) dw' 
Now in variables y 1 , ... , y n , w the equations (7.1) are written as 

(7.2i) = 

dy s V w 



Here we used the above notations (7.13) and (7.14). The equations (7.21) can be 
easily solved if we rewrite them as follows: 
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General solution of the equations (7.22) contains an arbitrary function of one vari- 
able h(w). It is given by the following formula: 

(7.23) a = h(w) V w {y\... ,y n ,w). 

Upon coming back to initial variables x 1 , ... , x n , v from (7.23) we obtain 

a = h(W(x 1 ,... ,x n ,v))x 

(7.24) 

x V w (x\... ,x n ,W(x\... ,x n ,v)). 

With the aim of additional verification we substitute the above expressions (7.24) 
and (7.15) into the equations (7.1). Let's do the appropriate calculations: 

^ = (h'(W) V w + h(W) V ww ) |^ + h(W) v sw , 

b s ^ = V s (h'(W) V w + h(W) V ww ) 

We add two above equalities and take into account formulas (7.16) and (7.17) for 
partial derivatives dW/dx s and dW/dv: 

(7.25) (J_ + hs ^ a = h[W )V sw . 

Similar calculations for the right hand side of the equations (7.1) yield 

(7.26) [a—\b s = h{W)V w V sw — . 

Comparing (7.25) with (7.26) and taking into account formulas (7.16) for the de- 
rivative dW/dv completes the proof of the following theorem. 

Theorem 7.2. Functions b\, . . . , b n and a satisfy nonlinear differential equations 
(7.1) and (7.2) if and only if they are determined by formulas (7.15) and (7.24). 

8. General formula for force field. 

Analyzing reduced normality equations above, we have found their general solu- 
tion. Now we are able to write formula for the force field of arbitrary Newtonian 
dynamical system admitting the normal shift on Riemannian manifold of the 
dimension n ^ 3. Let's substitute (5.5) into the formula for the components of 
force field F. Thereby we take into account that fields a and b determining scalar 
parameter A are fiberwise spherically symmetric: 

n 

(8.1) F k = a 7V fc + |v|]T& 4 (2N l N k -5l). 

i=l 

Formula (8.1) completely determines the dependence of force field on velocity vector 
v = !)N. Taking into account (7.15) and (7.24) we get 

F k = h(W) V w (x\... ,x n ,W)N k + 

(8.2) 



\v\J2v t (x\... ,x n ,W) (2N*N k -5l). 



i=i 
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Here h(W) is an arbitrary function of one variable; through Vi and V w we de- 
noted derivatives (7.13), where functions V{x 1 ,... , x n ,W) and W(x 1 ,... , x n ,v) 
are bound with each other by the relationships (7.12). 

Formula (8.2) for the force field of dynamical system admitting the normal 
shift contains the arbitrariness determined by one function of (n + 1) variables 
V(x 1 ,... ,x n ,w). Arbitrariness determined by the function h(W) can be elimi- 
nated by means of gauge transformation that changes a but doesn't change b: 

V(x\... ,x n ,w) -^V(x\... ,x n ,p-\w)), 

(8.3) W(x\... ,x n ,v) -^p(W(x\... ,x n ,v)), 

h{w)^h{p-\w)) p'{p-\w)). 

Hence transformation (8.3) doesn't change components of force field F, though it 
changes parameters V, W, and h in (8.2). If h(w) ^ 0, we can choose function 
p{w) such that h(w) p'(w) = 1. Upon doing gauge transformation (8.3) in this case 
we obtain the following formula for force field F: 

F k = V w (x\... ,x n ,W) N k + 

+ M^V^ 1 ,... ,x n ,W) (2N*N k -Sl). 

i=l 

Formula (8.4) is almost as universal as formula (8.2). The only exception is the 
case h = 0, which is not embraced by formula (8.4). 

9. Effectivization of general formula. 

Formulas (8.2) and (8.4) determine components of the force field of Newtonian 
dynamical system admitting the normal shift on Ricmannian manifold M. However, 
both these formulas have common fault. They are ineffective since we are to use 
function W(x 1 , . . . , x n ,v) determined implicitly by the equations (7.12). With the 
aim to get more effective formula we use the relationships (7.16) and (7.17). Let's 
rewrite these relationships as follows: 

(9.1) V w {x\... ,x n ,W(x\... ,x n ,v)) 



dW/dv ' 
dW/dx k 



(9.2) V k (x\...,x n ,W(x 1 ,...,x n ,v))- . )]v/()r . 

Substituting (9.2) into the formula (7.15), for the components of covector b we get 

dW(x\... ,x n ,v)/dx k 



(9.3) b k 



dWix 1 ,... ,x n ,v)/dv ■ 



Here W{x 1 , ... , x n ,v) can be understood as absolutely arbitrary function provided 
the derivative in denominator of the fraction in (9.3) is non-zero. The same function, 
upon substituting (9.1) into (7.24), determines the field a: 

„_ h(W(x\... ,x n ,v)) 



dW(x\... ,x n ,v)/dv' 
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Let's substitute the expressions (9.3) and (9.4) into the formula (8.1) for F: 

h(W(x\... ,x n ,v)) 
k dWix 1 ,... ,x n ,v)/dv k 

" dW(x\...,x^v)/dx* 
~ M ^ dW(x\...,x",v)/dv {2N Nk ~ h >- 

Let's rewrite formula (9.5) in terms of covariant derivatives. We formulate the 
result in form of the theorem. 

Theorem 9.1. Newtonian dynamical system on Riemannian manifold M of the 
dimension n ^ 3 admits the normal shift if and only if its force field F has the 
components determined by formula 

i— 1 

where W is fiberwise spherically symmetric scalar field from extended algebra of 
tensor fields on M with non-zero derivative 

i=l 

and h — h(W) is an arbitrary function of one variable. 

10. Kinematics of normal shift. 

Having explicit formula (6.6) for the force field of the dynamical admitting the 
normal shift, we are able to describe in details the process of normal shift of a 
given hypersurface S along trajectories of this dynamical system. According to the 
definition 3.2 one chooses some point po on S, then on some part S' — Os(po) of 
hypersurface S one should define the function v that determines the modulus of 
initial velocity on the trajectories of shift. At the point po this function v should 
be normalized by the condition 

(10.1) i/(po) = VQ, 

where v is some nonzero number (see condition (3.3) above). Let's choose local 
coordinates x 1 , . . . , x n on M in some neighborhood of the point p and local coor- 
dinates u , . . . , w" -1 on S in a neighborhood of the same point. The choice of local 
coordinates u 1 , ... , u 11 ^ 1 determines coordinate tangent vectors n, ... , r„_i to 
the hypersurface S. Let r l k be components of vector in coordinates x 1 , ... , x n 
on M. If hypersurface S is defined parametrically 

x 1 =x 1 (u\... ,u n ^), 

(10-2) 

x n =x n (u\... ,u n - 1 ), 
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then components of coordinate tangent vectors Ti, ... , r„_i are determined by 
derivatives of the functions x*(u 1 , . . . , w" _1 ) in (10.2): 



(10.3) 



du k ' 



Function v = v(p) = ^(u 1 , . . . ,u n x ) determines the initial velocity on the trajec- 
tories of shift starting from S (see formula (2.3) above): 

(10.4) v(t)| = i/(p)-n(p). 

lt=o 

Normal vector n(p) is determined up to a sign: n(p) — > ±n(p). Therefore without 
loss of generality we can assume that constant vq in (10.1) is positive. Then function 
v(p) in (10.4) is also positive. This means that 

(10.5) \v(t)\ =i/(p). 

t=o 

The condition (10.4) provides normality of shift at initial instant of time t = 
0. How to provide normality condition for other instants of time t ^ 0? For 
this purpose we consider the solution of Cauchy problem (2.1) for the system of 
differential equations (1.2), which describes Newtonian dynamical system with force 
field F. This solution is given by the functions (2.2). Let's write them as 

x 1 = x 1 (u 1 ,... ,u n ~ 1 ,t), 

(10-6) 

x n = x n (u 1 ,... ,u n -\t). 

Functions (10.6) describe the shift f t : S — ► St- In sufficiently small neighborhood 
S' = Os(po) of the point po and for sufficiently small values of t the map ft is 
a diffeomorphism: ft : S' — ► S' t . Due to this diffeomorphism we can carry local 
coordinates u 1 , . . . , from S to St- Then for any fixed t the functions (10.6) 

can be treated as parametric equations of hypcrsurfacc St similar to the equations 
(10.2). Let's calculate derivatives (10.3) for the functions (10.6). Doing this, we 
define vectors n, ... , r n _i tangent to all hypersurfaces S t : 

(10.7) T k = r k {u\... ,u n ~\ t). 

Denote by ipi the scalar product of the vector (10.7) and the vector of velocity: 

(10.8) Wfc = (r fe |v). 

In thesis [17] the functions tpi, . . . , ip n -i determined by formula (10.8) were called 
functions of deviation. Such functions play an important role in deriving the 
normality equations (3.1) and (3.2), since the condition of normality for the shift 
ft'- S — > St consists in identical vanishing of all functions (10.8): 

(10.9) <f k =<p k (u\... ,u n -\t) = 0. 
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From identical vanishing of the functions of deviation (10.9) it follows that their 
time derivatives at the initial instant of time t = are zero: 

(10.10) v» fc | =0. 

lt=o 

Moreover from the relationships (10.9) it follows that ipk for t — are zero as well: 

(10.11) <fk\ =0. 

lt=o 

In papers [6] and [7] it was shown that for Newtonian dynamical systems admitting 
the normal shift (in the sense of definition 3.2 and theorem 3.1) the conditions 
(10.10) and (10.11) are not only necessary, but also sufficient for identical vanishing 
of all functions of deviation (10.8). 

The conditions (10.11) are trivial consequences of (10.4). They give no informa- 
tion on how to choose the function v(p) on S. Therefore let's consider the conditions 
(10.10). Let's calculate ip by differentiating (10.8): 

(10.12) cp k = V t <fk = (V t r fe | v) + (x fe | F). 

In deriving (10.12) we took into account the equations of dynamics (1.2) written as 
V 4 v = F. Now let's calculate the derivative V t Tfc. It's the vector with components 



Q i ™ ™ 

(io.i3) v *^ = ^+EE r '. ' 

r=l s=l 



If we take into account formulas (10.3) determining r^, then from (10.13) we obtain 
(10-14) V ^ = £ + EE v r ^ = V„„'. 

r—1 s — 1 

In vectorial form (10.14) is written as V t Tfe = V u fcV. That is V t Tfc coincides with 
covariant derivative of vector function v(u 1 , . . . , t) with respect to parameter 
u h along fc-th coordinate line on hypersurface S t . Let's substitute the obtained 
expression for V*Tfe into the formula (10.12). This yields 

(10.15) 99 fe = (V„ fe v|v) + (r fe |F). 

Further, we take into account the obvious relationship V„k(v|v) = 2(V u fcv|v) 
and formula (10.5), which determines modulus of velocity vector for t = 0. Then 
we can bring formula (10.15) for the derivative ipk to the following form: 

dv 

(10.16) <p k \ =v^ + {r k \¥). 



\t=o 



Now, due to (10.16), the relationships (10.10) for derivatives (pk are written as 
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partial differential equations for the function v — v(u , . . . , u n ) on S: 
(10.17) ^ = -^(F\r k ). 

Let's substitute force field (9.6) into the equations (10.17) and take into account 
the fact that vector N for t = coincides with unitary normal vector n(p) on S. 
Upon rather simple calculations this yields 

(m8) ^ = -£w^- 

Let's multiply the equation (10.18) by W v and transfer the sum from left to right 
hand side of this equation. Moreover, let's write explicitly all derivatives: 

(10 19) — — V QW - 

dv du k dx % du k 

i=i 

It's not difficult to see that left hand side of (10.19) is the derivative of the func- 
tion W(x 1 , . . . , x n , v) with respect to u k upon substituting the functions (10.2) for 
x 1 , . . . , x n and the function u{v} , . . . , u™ _1 ) for v. Therefore the equations (10.19) 
are easily integrated in form of functional equation 

(10.20) W{x 1 {p) 1 ... ,x n (p),u(p)) = W Q = const, 

which determines the function v = v{p) = ^(m 1 , . . . ,u™ _1 ) in implicit form. The 
value of constant Wo in (10.20) is fixed by normalizing condition (10.1): 

W = W(x 1 (po),... ,x n (p ),is ). 

Theorem 10.1. In order to construct the normal shift of hypersurface S, given 
in parametric form by functions x 1 {p), ... , x n (p) from (10.2), along trajectories of 
Newtonian dynamical system with force field (9.6) one should determine the function 
u(p) in (10.4) by means of the equation (10.20). 

Having constructed the normal shift f t : S —> St along trajectories of dynamical 
system with force field (9.6), we obtain the family of hypersurfaces St. By changing 
the initial instant for counting the time t — > t + to we can treat each hypersurface of 
this family as initial hypersurface. Therefore on each of them the following equality 
similar to (10.20) is fulfilled: 

(10.21) W(p, |v|) = W (t) = const . 

Note that the values of constants Wo(t) in (10.21) can be different on different 
hypersurfaces St- Let's calculate the dynamics of Wo(t) in t. First find the time 
dynamics of the modulus of velocity vector on the trajectories of shift: 

(10.22) = Vt „ = = (N | F) = ±N k F k . 
dt v ^ 

k=i 
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Then substitute (9.6) into (10.22). As a result of this substitution we obtain 

dv _h(W) _^VtW dx* 

(W 6) dt~ w v ^ w v dt • 

i—l 

Now let's multiply (10.23) by W v and transfer the sum from left to right hand side 
of this equation. Moreover, let's write explicitly all derivatives: 

(1 °- 24) ^WW + W^ = h ^ 

In left hand side of (10.24) we see the time derivative of the scalar field on the 
trajectories of normal shift. Hence (10.24) is the required equation that determines 
time dynamics of constants Wo(t) in (10.21). Let's write this equation as follows: 



(10.25) = h(W ). 

If function h{w) in (9.6) is zero, then due to (10.25) the field W in (10.21) not 
only is constant on each separate hypersurface St , but has equal values on all such 
hypersurfaces. For h{w) ^ differential equation (10.25) is easily integrated. So, 
knowing the value of W on S, we can find its value on any one of hypersurfaces S t - 



11. Coordinates associated with Newtonian normal shift. 

It is known that the construction of geodesic normal shift of hypersurface S 
in Riemannian manifold M provides some special choice of local coordinates in 
a neighborhood of S. They are called semigeodesic coordinates (see [20] or 
[53]). Newtonian normal shift of hypersurface S also can provide some special 
choice of local coordinates in M. Let f t : S — > St be the normal shift of S along 
trajectories of Newtonian dynamical system with force field (9.6), function v{p) for 
which is fixed by normalizing condition (10.1) at the point po e S (without loss 
of generality we can assume that vq > 0). Then some neighborhood of the point 
po foliates into the union of not intersecting parts of hypersurfaces St- Choosing 
local coordinates u 1 , ... , u™" 1 on S in a neighborhood of po, we can carry them 
from S to St by means of shift diffcomorphism ft : S —* St- Therefore the set of 
n quantities u , ... , and t can be considered as local coordinates in M in 

some neighborhood of the point p . Such coordinates are called associated with 
normal shift f t : S — > S t - If we denote associated coordinates by a; 1 , ... , x n , then 
the functions x l (v}, . . . ,u n_1 ,t) in (10.6) become extremely simple: 

^(t* 1 ,... ,u n_1 ,t) = u\ 

(11.1) 

a;" -1 ^ 1 ,... ,u n -\t)=u n -\ 
x n {u\... ,u n -\t)=t. 



Differentiating functions (11.1) according to (10.3), we obtain the components of 
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vectors n, ... , r„_i tangent to the hypersurfaces S t in associated coordinates: 

{1 for i = k, 
for i^k. 

Differentiating functions (11.1) in t, we find components of velocity vector v: 

f 1 for i = n, 

"■2 v * = S\= ' 

^ tor i/n. 

Denote by ^ the modulus of velocity vector on trajectories of shift: 
(11.3) i/ = i/(u\... ,u n-1 ,£) = |v|. 

Doing this, we extend the domain of function v from (10.5), which is initially defined 
only for t = on initial hypersurface S. From (11.2) we obtain the relationship 



(11.4) 



V=\v\ = ^gnniu 1 , . . . ,U n ~ 1 ,t). 



On the other hand, the modulus of velocity vector can be calculated from the 
functional equation (10.21). One can solve this equation in explicit form by using 
the function V(x 1 , ... ,x n ,w) from (7.10): 



(11.5) 



|v| =V(u\... ,u n -\t,W Q (t)). 



Comparing (11.4) and (11.5), we get the following formula for diagonal component 
g nn of metric tensor in local coordinates associated with normal shift f t : S — > S t : 

(11.6) g nn = V(x\... ,x n ,W (x n )f- 

Following non-diagonal components of gij are zero due to normality of shift: 

(11.7) g nk = for all k=l, ... ,n - 1. 

Function of one variable Wo(t) from (11.5) and (11.6) is determined as the solution 
of ordinary differential equation (10.25) fixed by initial condition 



(11.8) 



W (t)\ 



W(p , vo). 



Initial condition (11.8) follows from (10.21), from the relationship (10.5), and from 
normalizing condition (10.1). 

The relationships (11.6) and (11.7) mean that matrix formed by components of 
metric tensor in local coordinates x 1 , ... , x n associated with Newtonian normal 
shift f t : S — ► S t has blockwise-diagonal structure: 



(11.9) 



G 



9 = 



g n 
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Trajectories of normal shift f t : S —> St correspond to the variation of parameter t 
in (11.1) by fixed values of parameters u 1 , . . . , u"^ 1 . Let's write the equation of 
Newtonian dynamics of points V t v = F for such trajectories. In associated local 
coordinates this vectorial equation reduces to the series of scalar equations. Taking 
into account the relationship (11.2), we obtain 

(11.10) lt„ = F fe , fc = l,...,n. 

Let's lower the upper index k in (11.10). Thereby we take into account blockwise- 
diagonal structure of matrix of metric tensor (11.9) and explicit formula (6.14) for 
components of metric connection. This yields 

(11.11) ^ = -2F k for k<n, 

(11.12) ^ = 2F n for k = n. 

The quantity g nn in left hand sides of (11.11) and (11.12) is determined by formula 
(11.6). Covariant components of force vector are determined by formula (9.6), or 
equivalent formula (8.2). By substituting (11.6) and (8.2) into (11.11) and into 
(11.12) we take into account that components of velocity vector on trajectories of 
shift are determined by formulas (11.2), while its modulus is determined by formula 
(11.3). Therefore for components of unitary vector N we have 

N k = jf y N k = \v\8%. 

If we remember all circumstances listed above, then by substituting (11.6) and (8.2) 
into the equalities (11.11) and (11.12) we find that these equalities turn to identities. 
So we get no restrictions for the choice of functions h(w) and V(x , . . . , x n , w). 
This is not surprising, since all restrictions due to normality of shift ft:S—*St 
are already handled by normality equations (3.1) and (3.2), and by explicit formula 
(8.2) that follows from these equations. As for the normalizing condition (3.3) for 
v, it is provided by the relationship (11.6), by the equation (10.25) for the function 
Wo(t), and by initial condition (11.8). 

The above result gives the answer to one of the questions by A. V. Bolsinov 
and A. T. Fomenko. It is formulated as follows: to what extent the normal shift 
of some particular hypersurface f t : S — > St characterizes the structure of force 
field of dynamical system admitting the normal shift? The answer is: yes, it 
characterizes, but partially; it doesn't determine F completely. 

Indeed, if normal shift ft : S — ► S t is already constructed, then constructing 
, . . . , x n in a neighborhood of S reduces to the choice of 
on initial hypersurface S. Diagonal component g nn 
of metric tensor in these coordinates determines the function of n variables 



(11.13) V(x\ ...,x n , W (x n )) = vW* 1 , • ■ • > * n ) 



(see the relationship (11.6) above). But by function (11.3) one cannot reconstruct 
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the function of n + 1 variables V(x 1 , ... , x n , w), which is contained in the formula 
(8.2) for components of force field F. 

If we suppose that h(w) — in formula (8.2), then from (10.25) it follows that 
Wo(x n ) — Wo — const. In this case formula (11.13) simplifies to 

(11.14) V{x\... ,x n ,W ) = y/g nn ( x i,... ,x n ). 

But for the fixed normal shift f t : S — > St the constant Wo is strictly fixed. Therefore 
by (11.14) we cannot reconstruct the function V(x 1 , . . . ,x n ,w) in whole. 

Now suppose that h(w) ^ 0. Let's consider gauge transformations (8.3) that do 
not change force field (8.2). Transformations (8.3) are supplemented by the rule 
for transforming Wo(t). It looks like 

(11.15) W (t) — p(W (t)). 

Gauge transformations (8.3) supplemented by the additional rule (11.15) preserve 
not only the force field F of dynamical system, but the maps of normal shift f t as 
well. Therefore they do not change the choice of associated local coordinates and 
the function g nn in right hand side of (11.13). Invariance of the left hand side of 

(11.13) with respect to these transformations is easily checked by direct calculations. 
At the expense of the gauge transformations (8.3) the case when h{w) ^ can be 
reduced to the case h(w) = 1 (see comment preceding formula (8.4)). For h(w) = 1 
by means of integrating (10.25) we get Wo(t) — Wo + 1, where Wo — const. This 
reduces (11.13) to the following form: 

(11.16) V{x\... ,x n ,W + x n ) = ^/g nn (x\... ,x«). 

Constant Wo in (11.6) is strictly fixed for the fixed normal shift f t :S—> S t - 
Therefore by (11.16) one cannot reconstruct the function (11.16) in whole. 

Consider a simple example. Let M = M™ be cuclidcan space with standard 
metric, and let S be hyperplane given by the equation x n = 0. Then the rela- 
tionships (11.1) determine parametrization of S and define the normal shift of this 
hypersurface f t : S — > S t , being the parallel displacement of S along n-th coordinate 
axis. Thereby g nn = 1. 

1. The above shift can be implemented by Newtonian dynamical system with 
identically zero force field F(r,v) = 0. This corresponds to the choice h(w) = 
and V(x 1 , . . . ,x n ,w) — w in formula (8.2), and to the choice Wo = 1 in formula 

(11.14) respectively. 

2. Function V(x 1 , . . . ,x n ,w) can be changed, keeping h(w) = and Wo = 1 
meanwhile. Let's set V(x 1 7 . . . , x n , w) = w + (1 — w) ■ fix 1 , . . . , x n ). Then the 
above normal shift of hyperplane S will be implemented by Newtonian dynamical 
system, force field of which is non-zero: 

F( r v) = 1 ~ |v| . 2(v|Vy)-v-|v| 2 -V^ 
1 — ip Ivl 



Here <p = <p(r) is an arbitrary function of coordinates x 1 , ... , x n , for which is 
natural to assume, that its values are distinct from 1. 
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3. Taking h(w) = 1, we can choose W(x 1 , . . . ,x n ) = w — x n in formula (8.2), 
and Wo = 1 in formula (11.16) respectively. For the force field F this yields 

„, . v 2 (v I M) • v - Ivl 2 • M 



The example, which we have just examined, confirms our conclusion that knowing 
the normal shift ft- S — > St of some particular hypersurface is not sufficient for to 
determine the force field of dynamical system implementing this shift, even if we 
know that this system belong to the class of systems admitting the normal shift. 

12. Blowing up the points. Generalization of theory 
as suggested by a. v. bolsinov and a. t. fomenko. 

Before now, studying normal shift ft : S — > St , we restricted ourselves to the case 
of smooth hypersurfaces and took parameter t small enough for hypersurfaces S t 
to be non-singular as well. However, one case with singularity appears to be inter- 
esting now. This is the case when hypersurface St collapses into a point at a time 
for some t = t . By reverting the direction of time we can speak about blowing up 
the point. Moreover, without loss of generality we can assume that to = 0. In this 
case we have singular initial hypersurface S = {po} consisting of only one point po, 

and a fan-shaped pencil of trajectories coming 
out from this point (see figure 12.1). Velocity 
vectors v = v(0) on these trajectories corre- 
sponding to the time instant t = belong to 
the tangent space T po (M). They determine a 
hypersurface s in the fiber of tangent bundle 
TM over the point p € M. It can be under- 
stood as "limiting variety" for hypersurfaces 
St in "infinitesimal scale" : 

(12.1) s = limA 

If initial values of velocity vectors on all tra- 
jectories at the point po are non-zero, then s 
have topology of (n — l)-dimensional sphere. The same topology is inherited by 
all hypersurfaces S t for sufficiently small values of parameter t. Let a be the unit 
sphere in the fiber of tangent bundle over the point po, let q be a point of this 
sphere, and let n(q) be radius- vector of the point q in T po (M). Then radius- vectors 
of the points on the hypersurface s are given by formula v(q) = v(q) ■ n(q), where 
v = v(q) is some positive function on unit sphere a, while trajectories coming out 
from the point po are determined by initial data 

(12.2) x k \ =x k (p ), x k \ = V {q)-n k {q) 

\t=0 \t=0 

for the equations of Newtonian dynamics (1.2). Components of force vector in (1.2) 
are determined by formulas (8.2). Due to normality of shift f t : S — > St all results 
obtained above in sections 10 and 11 remain valid for non-singular hypersurfaces 
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St with t = 0. Formula (11.5) from section 11 now is written as follows: 

(12.3) \v\ = V(x 1 (t,q),... ,x n (t,q),W (t)). 

We can return to initial form of formula (11.5) if we denote by u 1 , ... , u™ -1 local 
coordinates of the point q on unit sphere a. 

Modulus of initial velocity |v| = v(q), which is contained in formula (12.2), can 
be found by passing to the limit t — > in formula (12.3): 

(12.4) u(q) = V(x 1 ( Po ),... ,x n ( Po ),Wo). 

Here x (po), • ■ ■ , x n (po) are coordinates of the point po, and Wo is initial value of 
function Wo(t) for t = 0. The function Wo(t) itself is determined as solution of 
ordinary differential equation (10.25). Note that right hand side of (12.4) doesn't 
depend on q, i. e. function v{q) is constant: 

(12.5) v(q) = V( Po , W ) = vo = const . 
Constant Wo can be expressed through v if we take into account (7.12): 

(12.6) Wo = W(po, vo) = const . 

On account of the relationships (12.5) and (12.6) we can rewrite (12.2) as follows: 

(12.6) x k \ =x k (po), x k \ = v -n k (q). 

lt=o lt=o 

Thus, we can state a theorem that follows from the results of sections 10 and 11 by 
passing to the limit (12.1). 

Theorem 12.1. Suppose that on Riemannian manifold M some Newtonian dy- 
namical system admitting the normal shift is defined, i. e. we have the system with 
force field F given by formula (8.2). Then for any point po £ M and for any positive 
constant vo > initial data (12.6) determine the normal blow-up f t : po — > St of 
the point p along trajectories of this dynamical system. 

Consideration of normal blow-ups for separate point of the manifold M gives the 
opportunity for further development of the theory of dynamical systems admitting 
the normal shift. We can formulate the following definition similar to definition 3.1. 

Definition 12.1. Newtonian dynamical system (1.2) with force field F on Rie- 
mannian manifold M is called a system admitting the normal blow-ups of 

points if for any point po G M and for arbitrary positive constant vq > ini- 
tial data (12.6) determine the normal blow-up f t : p — > S t of the point p along 
trajectories of this dynamical system. 

The idea of constructing new theory on the base of definition 12.1 was suggested 
by A. V. Bolsinov and A. T. Fomcnko in February of 2000 in the seminar at Moscow 
State University during the discussion on the results of thesis [17]. Theorem 12.1 
shows that dynamical systems with force field (8.2) are included into the framework 
of new theory. But, possibly, one can find some new dynamical systems, which 
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aren't embraced by formula (8.2). We can compare this situation with that of the 
theory of distributions, where narrowing class of test functions extends the class 
of distributions. Here, narrowing class of initial hypersurfaces in the construction 
of normal shift f t : S — > St to singular one-point sets, we have a good chance to 
extend class of dynamical systems that can implement such shift. Is it really so? 
The answer to this question can be given only as a result of constructing new theory. 
But this falls out of the limits of this paper. 

13. On the problems of metrizability. 

Problem of metrizability arose on initial stage of developing theory of dynamical 
systems admitting the normal shift by testing theory for non-triviality. It was noted 
(see papers [5] and [7]) that if metric g is conformally equivalent to the basic metric 
g of Riemannian manifold M, i. e. if we have 

(13.1) g = e- 2 'g, 

then geodesic flow of metric g is a dynamical system admitting the normal shift 
with respect to metric g. Its force field is given by formula 

n 

(13.2) F k = -\v\ 2 V k .f + 2Y,Vsfv s v k . 

Here V/ is a gradient of scalar field / determining conformal factor e~ 2 ^ in (13.1). 

Definition 13.1. Newtonian dynamical system on Riemannian manifold is called 
metrizable system if it inherits trajectories of the system with force field (13.2). 

Trajectory inheriting and trajectory equivalence for two dynamical systems are 
understood in the sense of the following definitions. 

Definition 13.2. Suppose that on the Riemannian manifold M two Newtonian 
dynamical systems are defined with force fields F and F respectively . Say that 
second system inherits trajectories of the first system if any trajectory of the 
second system as a line (up to a regular reparametrization) coincides with some 
trajectory of the first system. 

Definition 13.3. Two Newtonian dynamical systems on the Riemannian manifold 
M are called trajectory equivalent if they inherit trajectories of each other. 

Note that definitions 13.2 and 13.3 are somewhat different from corresponding 
definitions used in papers [54-60] . Our definitions are more specialized and adopted 
to the case of Newtonian dynamical systems with common configuration space. 

Metrizable dynamical systems are trivial regarding to their use in the construc- 
tion of normal shift. Normal shift along trajectories of such systems, in essential, is 
reduced to geodesic normal shift. Therefore in papers [5] and [9] we considered the 
problem of describing all metrizable Newtonian dynamical systems admitting the 
normal shift, and in paper [15] we constructed examples of non-metrizable ones. 
Main result of papers [5] and [9] is formulated in the following theorem. 
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Theorem 13.1. Newtonian dynamical system admitting the normal shift on the 
Riemannian manifold M with metric g is metrizable by means of conformally equiv- 
alent metric g = e~ 2 f g if and only if its force field is given by formula 



where H = H(v) is some arbitrary function of one variable. 

Theorem 13.1 solved the problem of describing dynamical systems admitting the 
normal shift and being metrizable by means of conformally equivalent metric. How- 
ever, the requirement of conformal equivalence of metrics g and g in this theorem 
is a priori. One can exclude this requirement. Then geodesic flow of metric g will 
correspond to the dynamical system with less special force field 



in the metric g. One can step further, i. e. one can avoid metric g at all assuming 
Iy- in (13.3) to be components of some symmetric affine connection in M. 

Definition 13.4. Newtonian dynamical system on Riemannian manifold M is 
called metrizable by geodesic flow of affine connection T in M if it inherits 
trajectories of this geodesic flow. 

This is the very treatment of the concept of metrizability that was considered 
in paper [11], which is not published unfortunately. In paper [11] the following 
theorem was proved. 

Theorem 13.2. Newtonian dynamical system admitting the normal shift on the 
Riemannian manifold M of the dimension n ^ 3 is metrizable by geodesic flow 
of affine connection T if and only if this connection is (at least locally) a metric 
connection for some metric g = e~ 2 ^ g, which is conformally equivalent to basic 
metric g of the manifold M. 

In other words, if there exists some Newtonian dynamical system which inherits 
trajectories of geodesic flow of affine connection T, and which is admitting the 
normal shift, then connection T is necessarily a metric connection for some metric 
g = e~ 2 ^ g. And conversely, if T is defined by metric g = e -2 ^ g, then one can 
find some dynamical system admitting the normal shift and inheriting trajectories 
of geodesic flow of T. Though the converse proposition of the theorem is obvious, 
since geodesic flow of metric g = e~ 2 ^ g admits the normal shift in metric g. 

Theorem 13.2 shows that a priori assumption on conformal equivalence of metrics 
g and g used in papers [5] and [9] at the first approach to the problem of metriz- 
ability doesn't cause the loss of generality. Proof of the theorem 13.2 is based on 
the following fact, which was proved in [5]. 

Theorem 13.3. Suppose that force field F of the first Newtonian dynamical system 
on Riemannian manifold M is a homogeneous function of degree 2 with respect to 
components of velocity vector v in the fibers of tangent bundle TM. Then second 



n 



H{ve f)ef 

M 



F fe = -|v| 2 V fe / + 2^V s /« s v k + 



n 



n 



(13.3) 




32 



RUSLAN A. SHARIPOV 



Newtonian dynamical system inherits trajectories of the first system if and only if 
its force field F is given by the following formula: 

F(p,v)=F(p,v) + fe^.v. 

Components of force field (13.3) are quadratic functions with respect to the com- 
ponents of velocity vector. Hence the for force field of dynamical system inheriting 
trajectories of geodesic flow of affinc connection T we have the formula 

n n 

(13.4) F k = E E^ - r-jK v> + H N k . 

i=i j=i 

Further proof of theorem 13.2 in unpublished paper [11] consisted in substituting 

(13.4) into the normality equations (3.1) and (3.2). Here we shall give more simple 
proof of this theorem based on comparison of formulas (13.4) and (8.2). Let's denote 

= T k j — T k y The quantities M k j are the components of some (not extended) 
tensor field M on the manifold M. It is called the field of deformation or the field 
of variation for connection L. Now formula (13.4) is written as follows: 

n n 

(13.5) F k = -J2J2 M ii vivj + HNk - 

»=i j=i 

Proof of the theorem 13.2. Let's contract both sides of formula (13.5) with com- 
ponents of orthogonal projector P from (4.3). This immediately excludes the term 
containing scalar function H = H(x 1 , ... , x n , v 1 , . . . , v n ): 



n n n 



(i3.6) E p fe Ffe = -EEE n m% v* vk 

fe=l fe=l i=l j=l 

Similar contracting in formula (8.2) cancels the entry of function h(W): 

n n 

(is.?) e p i pk - -m 2 E p l uk{ - x ^ ■■■> xn > w )- 

fe=l fe=l 
Here in (13.7) by U k — U k (x 1 7 . . . ,x n , |v|) we denote the following quantities: 

(13g) uk _£ g ki V i (x\... ,x»,W(x\... ,x n M)) _ 

i=l 

From (3.10) we sec that U k are components of extended vector field U, which 
is fiberwise spherically symmetric (see definition 5.1 above). Let's compare the 
relationships (13.7) and (13.6). Their left hand sides coincide. Hence we can equate 
right hand sides of these two relationships: 

n n n n 

( 13 - 9 ) E E E n M k v* v> = | V | 2 e n u k . 

fc=i i=i j=i fe=i 
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This relationship is remarkable, since the dependence on v in it is almost explicit. 
Indeed, the quantities do not depend on v, while quantities U k depend only on 
modulus of velocity vector. Let's fix the point p <G M. This means that we fix local 
coordinates x 1 , ... , x n . Then rewrite the relationship (13.9) in vectorial form: 

(13.10) P(M(v, v) - |v| 2 • U) = 0. 

Here M(v,v) is vector valued quadratic form determined by tensor M when twice 
contracting it with vector v. Note that according to (13.10) the operator of projec- 
tion P, when applied to the expression M(v,v) — |v| 2 • U, yields zero. Therefore 
the equality (13.10) can be rewritten as follows: 

(13.11) M(v,v) - |v| 2 -U(|v|) = A(v) - v. 

Vector U = U(|v|) in (13.11) depend only on modulus of velocity vector, while 
scalar A = A(v) can contain full scale dependence on v. Let's study this dependence. 
Consider vectors U(|v|) and U(|a • v|), where a is a number. Remember that in 
theorem 13.2 we deal with multidimensional case n ^ 3. In the space of the 
dimension n > 3 vector v ran be turned so that it doesn't belong to the linear span 
of vectors U(|v|) and U(|a • v|), while its modulus |v| being preserved unchanged. 
Let's substitute a ■ v for vector v into the equality (13.11): 

(13.12) a 2 ■ M(v, v) - a 2 |v| 2 • U(|a • v|) = a \(a ■ v) • v. 

Then multiply both sides of (13.11) by a 2 and subtract the obtained equality from 

(13.12) . As a result we get the following relationship: 

(13.13) a|v| 2 • (U(|a • v|) - U(|v|)) + (A(a • v) -aA(v)) • v = 0. 

For the case when vector v doesn't belong to linear span of vectors U(|v|) and 
U(|a • v|) from the equality (13.13) we derive 

(13.14) U(|a-v|)=U(|v|), 

(13.15) A(a-v) =qA(v). 

Though the equality (13.14) holds for the case when v belongs to linear span of 
U(|v|) and U(|a • v|) as well, since U depend on |v|, but not on the direction of 
v. Substituting (3.14) back to (3.13), we prove (3.15) for all v ^ 0. For v = the 
value of A(v) is not determined by formula (13.11). Therefore we can extend the 
function A(v) by taking A(0) = 0. This cancels the restriction v ^ in applying 
formula (13.15). 

Due to (13.14) vector U do not depend on v at all. Therefore left hand side of 
(13.11) is quadratic function in v. The equality (13.11) can be rewritten as 

(13.16) K(v) = K(v, v) = A(v) • v. 

Any quadratic in v function satisfies the following identity, which can be checked 



3d 



RUSLAN A. SHARIPOV 



by direct calculations: K(vi +v 2 ) +K(vi — v 2 ) = 2-K(vi)+2-K(v 2 ). Substituting 

(13.16) into this identity, we get the equality for A(v): 

A(vi + v 2 ) • (vi + v 2 ) + A(vi - v 2 ) • (vi - v 2 ) = 2 A(vi) • Vi + 2 A(v 2 ) • v 2 . 

If vectors vi and v 2 are linearly independent, then the above vectorial equality 
leads to the pair of scalar equalities: 

A(vi + v 2 ) + A(vi - v 2 ) = 2 A(vi), 
A(vi + v 2 ) - A(vi - v 2 ) = 2 A(v 2 ). 

Let's add them and divide the result by 2. Then we get the relationship 

(13.17) A(vi +v 2 ) = A(vi) + A(v 2 ). 

If vectors vi and v 2 are linearly dependent, then (13.17) follows from (13.15). The 
relationships (13.15) and (13.17) mean that A(v) is a linear function in v. 

Thus A(v) = (v | A), where A is some (not extended) vector field on M. When 
applied to quadratic form M(v,v) in (13.11), this yields 

(13.18) M(v,v) = |v| 2 -U+(v|A)-v. 

Above we have proved that U doesn't depend on v. Therefore U is also some (not 
extended) vector field on M. Now let's return to (13.8) and rewrite this equality 
in terms of covariant components of the vector U: 

Ui(x\... ,x n )\v\=V i (x 1 ,... ,x n ,W{x 1 ,... ,x n ,\v\)). 

The quantity |v| plays the role of independent variable in this equality. Let's 
substitute V(x 1 , . . . ,x n ,w) for |v| and take into account the relationships (7.12): 

U^x 1 ,... ,x n )V(x\... ,x n ,w) = V i (x 1 ,... ,x n ,w). 

Now remember the relationships (7.13). They show that Uiix 1 , . . . ,x n ) is a loga- 
rithmic derivative of the function V(x 1 , . . . , x n , w) with respect to the variable x l . 
This logarithmic derivative doesn't depend on w: 

(13.19) ^ ^ ) = g]fe^) . 

From (13.19) it follows that there exist (at least locally) two functions, a function 
/ = f{x 1 , . . . , x n ) and a function p = p(w) such that 

df(x\...,x n ) 

(13.20) dx l 

V = cxp(f(x\... ,x n ))p(w). 
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At the expense of gauge transformation (8.3) function V = V(x 1 ,... ,x n ,w) in 

(13.20) can be brought to the following form: 

(13.21) V = exp(f(x 1 ,... ,x n ))w. 

Let's substitute (13.21) into the formula (8.2) for the components of force field F: 

n 

(13.22) F k = fc(|v| e-f) e? A, + ]T M 2 U> (2 N* N k - 5%) . 

i=i 

Now let's compare formula (13.22) with formula (13.5) and take into account the 
above relationship (13.18) for tensor field M. This yields 

(13.23) g ' M | v V | l 2 e ' /)e ^ (2U + A|N). 

Here N is unitary vector directed along the vector of velocity. Note that left hand 
side of (13.23) depend only on modulus of velocity vector v, while right hand side 
depend on the direction of this vector. Therefore from (13.23) we get 

(13.24) H = h(\v\e~ f )e f , A = -2U. 

The relationships (13.20) and (13.24) completely determine the components of ten- 
sor field M. For components of connection V we have 

n 

(13.25) = I* - Vif5] - VjfS* +Y^9 kq V,/ //,... 

Substituting (13.25) into (13.3), we come to the force field (13.2). Force field (13.2) 
corresponds to geodesic flow of metric connection for metric g = e~ 2 f g. Thus, 
theorem 13.2 is completely proved. □ 
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